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SUMMARY 

NASA Lewis Research Center Is currently developing probabilistic struc- 
tural analysis methodology for select Space Shuttle Main Engine (SSME) com- 
ponents. This methodology consists of the following program elements: 

(1) composite load spectra, (2) probabilistic structural analysis methods, 

(3) probabilistic finite element theory - new variational principles, and 

(4) probabilistic structural analysis application. The methodology has led to 
significant technical progress In several Important aspects of probabilistic 
structural analysis. The program and significant accomplishments to date are 
summarized In this paper. 


INTRODUCTION 

It Is becoming Increasingly evident that deterministic structural analysis 
methods will not be sufficient to properly design critical structural compo- 
nents for upgraded Space Shuttle Main Engines (SSME). Structural components 
In the SSME are subjected to a variety of complex, severe cyclic and transient 
loading conditions Including high temperatures and high temperature gradients. 
Most of these are quantifiable only as best engineering estimates. These com- 
plex loading conditions subject the material to coupled nonlinear behavior 
which depends on stress, temperature, and time. Coupled nonlinear material 
behavior Is nonuniform, Is very difficult to determine experimentally, and 
perhaps Impossible to describe deterministically. In addition critical SSME 
structural components are relatively small. Fabrication tolerances on these 
components, which In essence are small thickness variations, can have signifi- 
cant effects on the component structural response. Fabrication tolerances by 
their very nature are statistical. Furthermore the attachment of the compo- 
nents to the structural system generally differs by some Indeterminant degree 
from that which was assumed for designing the component. In summary, all four 
fundamental aspects - (1) loading conditions, (2) material behavior, (3) geo- 
metric configuration, and (4) supports - on which structural analyses are 
based, are of a statistical nature. One direct way to formally account for all 
these statistical aspects Is to develop probabilistic structural analysis 
methods where all participating variables are described by appropriate proba- 
bilistic functions. 

NASA Lewis Research Center Is currently developing probabilistic struc- 
tural analysis methods for select SSME structural components. Briefly, the 
deterministic, three-dimensional. Inelastic analysis methodology developed 
under the Hot Section Technology ((HOST) and R&T Base programs) Is being 
augmented to accommodate the complex probabilistic loading spectra, the ther- 
movlscoplastlc material behavior, and the material degradation associated with 
the environment of space propulsion system structural components representative 



of the SSME, such as turbine blades, transfer duct, and liquid-oxygen posts 

(fig. 1). 

The development of probabilistic structural analysis methodology consists 
of the following program elements: (1) composite load spectra, (2) probabi- 

listic structural analysis methods, (3) probabilistic finite element theory - 
new variational principles, and (4) probabilistic structural analysis 
application. The development of the probabilistic structural analysis meth- 
odology Is a joint effort of NASA Lewis In-house research, contracts, and 
grants. The objective of this paper is to describe briefly the program and to 
highlight significant accomplishments to date. 


PROGRAM DESCRIPTION/PARTICIPANTS 

The major part of the program consists of two multiyear contracts. One 
contract (NASA contract NAS3-24382) Is for the development of composite load 
spectra. The prime contractor Is Rocketdyne, a division of Rockwell Interna- 
tional Corporation, with Battelle Columbus Laboratory, as a subcontractor. 
Rocketdyne Is responsible for the overall program and Integrated computer codes 
while Battelle Is responsible for developing probabilistic models. The other 
contract (NASA contract NAS3-24389) Is for the development of probabilistic 
structural analysis methods. The prime contractor Is Southwest Research 
Institute with several subcontractors. The Initial participants on this con- 
tract are: Southwest Research Institute - project management and solution 

strategies development; MARC Analysis Research Corporation - code development 
for probabilistic finite element methods; Rocketdyne Division, Rockwell Inter- 
national Corp. - Space Shuttle Main Engine (SSME) design and hardware exper- 
ience; Prof. Paul Wlrschlng (University of Arizona) - probabilistic and 
reliability methods; Prof. Gautam Dasgupta (Columbia University) - stochastic 
finite elements; and Prof. Satya Atlurl (Georgia Tech) - advanced material 
constitutive models. The remaining part of the program consists of one grant 
(NAG 3-535) with Northwestern University and the In-house effort which Is sup- 
ported by support service contract personnel from Sverdrup Techntflogy Incor- 
porated under NASA contract NAS3-24105. The program Is summarized In figure 2. 

The research activities have led to significant technical progress In 
several Important aspects of probabilistic structural analysis. Technical 
progress Is reported In monthly reports and oral contract progress reviews. 

In addition, progress was reported at the Structural Integrity and Durability 
of Reusable Space Propulsion Systems Conference held at NASA Lewis Research 
Center on June 4-5, 1985 (refs. 1 to 7). Furthermore, progress was reported 
In the ASME Winter Annual Meeting In Miami, Florida, November 1985, where a 
whole session was devoted to this subject (refs. 8 to 12). The information 
presented In this paper summarizes progresses reported In all these sources. 


COMPOSITE LOAD SPECTRA 

The main thrust of the composite load spectra (CLS) contractual effort Is 
to develop generic probabilistic models for the various Individual loads, their 
probable combinations for the select SSME components and attendant computer 
codes (refs. 2 and 3). The tasks of the contractural effort are summarized In 
figure 3. The type of available Information to aid In the development of the 
generic probabilistic models Is summarized In table 1 for the high pressure 
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fuel turbopumps (HPFTP) blade. Limited but typical measured data for compar- 
able Information Is shown In figure 4. The logic chart for analyzing and 
reducing test or flight measured data Is shown In figure 5. The types of 
Individual loads Identified to date for four SSME components are summarized In 
table 2 where the source for obtaining data for these Individual loads Is also 
shown. As can be seen In this table each component Is subjected to combina- 
tions at least nine Individual load conditions. 

The generic probabilistic models for each of the Individual loads consist 
of 4 parts: (1) a steady state load, (2) a periodic load, (3) a random load, 

and (4) a local spike. Each of these parts Is formally described with a mean 
and a variance or a standard deviation. The shape of each part Is determined 
using three different probabilistic distributions. Three probabilistic methods 
are also used for the phasing of these four parts Into a single probabilistic 
load model. It Is believed at this time that the four parts for each load 
condition will be sufficient to realistically represent present and anticipated 
Individual loading conditions. The shape of each of the four parts Is 
described parametrically with undetermined coefficients. These coefficients 
are selected by combinations of known data, computationally simulated data, and 
Information provided by experts. 

The generic model for the time phasing (or combination) of the Individual 
loads Into the composite load spectra Is also formulated using three probabi- 
listic methods. The formulation Is based on combinations of limited known/ 
anticipated data (of the type shown In fig. 4 for example). Information esti- 
mated by experts, and largely on the probabilistic synthesis of probabilist- 
ically occurring events. The rationale for using various levels of progressive 
sophistication of probabilistic models for both the Individual loads and the 
composite load spectra Is to "balance" the uncertainties associated with the 
various estimates used In the formulations of these models. The models are 
then valldated/adapted by using appropriate structural analyses, and assessing 
the resulting structural responses. The selection of the structural analyses 
to be used Is obtained from point design Information which Is either known or 
computationally simulated by using expert opinions. Structural analyses of 
these types are also used to determine the reliability and the respective level 
of confidence for the composite load spectra to be applied to a specific 
component. 

There are four Important considerations In the development of each model: 

(1) the ability of the model to handle nonstandard distributional forms, 

(2) the treatment of nonstationary processes, (3) the handling of physical 
dependencies In the model, and (4) the ability of the method to operate effi- 
ciently so that It will be able to be Included In an expert system computer 
code. The structure for each probabilistic generic model has a distribution 
fitting routine for the Individual loads with a barrier crossing method. Dis- 
crete Probability Distribution (DPD) method, and a Monte Carlo method for the 
combined and composite load models. A transformation method for changing non- 
stationary load processes Into stationary processes Is also Included In the 
model so that the barrier crossing techniques can be used for a broader spec- 
trum of problems. Because these models are based on the probability density 
function representing the frequency of the various load levels occurring, all 
forms of load shape curves can be handled Including nominal (rectangular 
pulse), periodic, periodic over nominal, and random over nominal. The case of 
spike (transient) loads occurring Is handled with a simulation method, or. If 


3 



appropriate, the barrier crossing method. A typical Individual load record 
determined using this approach Is shown In figure 6. 


EXPERT SYSTEM COMPUTER CODE 

The Integrated computer code for the composite load spectra Is an execu- 
tively driven modular software system that will Incorporate the various Indiv- 
idual and composite load spectra models. This code Is configured as an "expert 
system". Modules In this code will contain the various probabilistic models 
with data dependent coefficients and/or functions and with the Increasing 
levels of sophistication with respect to prediction reliability and confidence 
level. Associated with these models are guidelines to select the coefficients 
and functions In the generic models to achieve a specified prediction require- 
ment. The potential complexity and the expertise Inherent In generic models 
to construct specific load spectra, justifies embedding these models In an 
Integrated software system configured as an expert system that can advise users 
In the employment of the generic load models. 

The code Is being developed In four Incremental versions. Each subsequent 
version will add a new engine component, additional load types and more sophis- 
tication to the probabilistic load definition and decision making process. The 
technical core of the final Integrated system will be an Interactive expert 
system that will: (1) construct specific load spectra models based on user 

supplied descriptions of the component and the load environment; (2) Incorpo- 
rate the probabilistic models so as to enable a dialogue between the system and 
the user which will help the user select the best model parameters for his 
problem; and (3) be able to describe the process of constructing a specific 
load spectra model. Including decisions made by the expert system and the 
rationale for those decisions. The key features of the code are summarized In 
figure 7. 

The expert system Is being configured so that specific simulated load 
spectra models are built by accessing a knowledge base of facts and rules. 

The knowledge base contains a module of decision-making data (facts) and a 
module of rules and decision criteria (rules) for constructing the load spectra 
model. All numerical results needed to construct the specific load spectra 
model will be computed or retrieved automatically by the expert system. The 
knowledge base and the Inference mechanisms needed to search the knowledge base 
will be developed as part of this project. They will be strongly dependent 
upon the probabilistic generic load spectra models and engine component 
loading. A schematic of the structure of the expert system Is shown In 
figure 8. 

The expert system Is being constructed to easily accommodate additions and 
deletions to the knowledge base. This allows the expert system to "get 
smarter" and to adapt to new Information, both In the test and validation stage 
of each code version and In the expansion of the expert system to Incorporate 
future models In the later versions of the code. The expert system Is Inter- 
active to allow the system to query and guide the user as to the system opera- 
tion (e.g., required user's Input). The knowledge base and Inference mechanism 
Is being designed to minimize redundant data requests from the system to the 
user. The software system Is modular In conformance with modern programming 
practice. The source program Is written In the FORTRAN language In order to 
assure a stand alone and portable code. 
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PROBABILISTIC STRUCTURAL ANALYSIS METHODS (PSAM) 

The focus of the PSAM contractual effort (refs. 4, 5, and 9 to 21) Is to 
develop analysis methods and computer programs for predicting the probabilistic 
response of critical structural components for current and future space pro- 
pulsion systems. This methodology will play a central role In establishing 
Increased system performance and durability. The Initial activity of PSAM Is 
the development of a probabilistic finite element code for the probabilistic 
structural analysis of the select SSME components. The computer code being 
developed Is NESSUS (Numerical Evaluation of Stochastic Structures Under 
Stress). It Is based on the Integration of existing technologies (fig. 9). 

The development of the NESSUS code Is scheduled to take 3 years. First 
year efforts Involve the formulation of the probabilistic analysis strategy and 
the development of a probabilistic linear analysis code. The ultimate goal of 
the 3-year program Is the development of a finite element code capable of per- 
forming nonlinear dynamic analysis of structures having stochastic material 
properties, geometry, and boundary conditions and subjected to random loading 
(fig. 10). 

Three levels of sophistication are pursued for the stochastic description 
of the structural problem namely: 

Level 1: Homogeneous random variable for stiffness, mass, damping, and 
external loading 

Level 2: Stochastic characterization of variables at the element level, 
with specified Interelement correlations 

Level 3: Stochastic Interpolation of variables within a finite element 

Two alternative probabilistic analysis methods are being developed, 
allowing for all three levels of modeling sophistication. These two methods 
are: 


(1) A probability Integration method, providing a direct estimate of the 
reliability of the structural configuration under study. 

(2) A simulation method, providing the means to verify the results obtained 
with method 1. 

The finite element library of NESSUS Is shown In figure 11 together with 
the perturbation variables and the processors. A demonstration problem of the 
NESSUS status Is briefly described below to Illustrate PSAM-type of computa- 
tional simulation. The problem chosen Is a curved shell representing a HPFTP 
blade. The finite element model for this shell-blade consists of 48 bilinear 
shell elements, 63 nodes, and 378° of freedom. The shell Is loaded with random 
pressure and temperature fields. These fields have a constant mean value over 
the surface but are spatially correlated through an exponential decay function. 
For pressure, the correlation was assumed strong In the spanwlse and weak In 
the chordwlse directions. The opposite was assumed true for the temperature 
field. 

The other random variables were shell thickness and stiffness at the base. 
The thickness variability could come from manufacturing processes, while base 
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stiffness variations may arise from normal assembly procedures. Proper defin- 
ition of boundary conditions as reflected In the base stiffness may not be well 
defined and could be a major source of uncertainty In a structural analysis. 
This would be especially true for an eigenvalue analysis where natural fre- 
quency Is the critical response variable. 

Randomness was Included In the elastic modulus and coefficient of thermal 
expansion of the material by assuming they were deterministic functions of 
temperature. For the strength analysis, material yield strength was also 
assumed to be stochastic. The distributions, their means and coefficients of 
variation for the random variables, were not taken to represent a particular 
SSME material, structure, or operating environment. They were selected as 
being representative of what might be expected. The higher coefficients of 
variations for the base stiffness factor and the pressure reflect the greater 
uncertainty these random variables are expected to have In actual service. For 
this analysis the only limitation on the choice of distributions Is that the 
temperature and pressure fields were taken to be normal. 

The NESSUS code was executed for a number of preselected perturbations 
about the deterministic state. In this example, the deterministic solution was 
defined at the mean value of the random variables. A total of 40 perturbed 
solutions were obtained at +1.5 and +3.0 standard deviations of the 11 random 
variables required for the combined stress performance function (Von Mlses 
criterion). The use of the modified Newton-Raphson Iteration method allowed 
for efficient computation of the perturbed solutions without repeated reformu- 
lation and resolution of the matrix equations. Clearly the number of pertur- 
bations required to construct a performance function depend on the number of 
random variables; consequently, for computational efficiency, only the signi- 
ficant variables should be retained In the analysis. The results obtained for 
excedence of Von Mlses type stress are summarized In figure 12 (lower right) 
where the blade-shell finite element model, the NESSUS code flow chart and the 
Input random variables are also shown. 


VARIATIONAL THEORY FOR PROBABILISTIC FINITE ELEMENTS 

A small but Important part of the probabilistic structural analysis meth- 
odology Is the development of variational principles for formulating probabi- 
listic finite elements. This part Is considered to be fundamental and Is 
pursued under a grant (refs. 6 to 8) which focuses on embedding the probabi- 
listic aspects In a variational formulation. A variational approach to prob- 
abilistic finite elements enables It to be Incorporated within standard finite 
element methodologies. Therefore, once the procedures have been developed, 
they could easily be adapted to existing general purpose programs. The varia- 
tional basis for these methods enables them to be adapted to a wide variety of 
structural elements and to provide a consistent basis for Incorporating prob- 
abilistic features In many aspects of the structural problem: (1) displace- 

ments, (2) boundary conditions, (3) body forces resulting from acceleration 
loads and, (4) any other features that cannot be clearly established. For 
example, the well known dilemma as to whether a shell Is clamped or simply- 
supported at a boundary, could also be treated more rationally by using a 
probabilistic distribution for this boundary condition. Relevant progress of 
this research effort Is summarized below. 
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A methodology has been completed which can embed the probabilistic distri- 
bution of the constitutive properties and loads (l.e., material uncertainties 
and load uncertainties) with a finite element variational approach. The cor- 
responding probabilistic distribution of the elemental nodal forces Is 
assembled Into a description of the probabilistic distribution of the nodal 
forces for the complete model. The appropriate mean structural response due 
to these uncertainties Is then determined efficiently. The basic concept of 
this approach Is to Incorporate the probabilistic distributions, as reflected 
In the variance, of the material properties and the loading conditions to 
obtain the corresponding variances In the elemental nodal forces of the finite 
element model. On the basis of the variance In the elemental nodal forces, the 
variance In the final solution Is determined In the usual deterministic solu- 
tion procedures. 

Efficient numerical algorithms were developed for obtaining the probabi- 
listic sensitivity element matrices which reflect the effects of randomness on 
response variables such as displacements, stresses, etc. The randomness Is due 
to the preassigned probabilistic descriptions of the material properties and 
loads. 

A pilot computer code was completed. Solutions obtained using this code 
have been compared to the Monte Carlo methods and the Hermlte Gauss Quadrature 
Integration schemes. The cost of the new method Is substantially lower. This 
was demonstrated with a ten-bar probabilistic nonlinear system where the random 
variables are the yield stresses. Exploitation of this characteristic In any 
probabilistic analysis may offer significant savings and warrant further 
Investigation. A flow chart of the pilot computer code Is shown In figure 13 
together with some representative results. 


PROBABILISTIC STRUCTURAL ANALYSIS OF SSME TURBOPUMP BLADES 

A probabilistic study has been Initiated In-house at NASA Lewis Research 
Center (ref. 7). The first objective of this study Is to evaluate the 
geometric and material properties tolerances on the structural response of 
turbopump blades. During this study, a number of Important probabilistic var- 
iables have been Identified which are considered to affect the structural 
response of the blade. In addition, a methodology has been developed to sta- 
tistically quantify the Influence of these probabilistic variables In an 
optimized way. The Identified variables Include random geometric and material 
properties perturbations, different loadings and a probabilistic combination 
of these loadings. Influences of these probabilistic variables are quantified 
by evaluating the blade structural response. The structural response variables 
Include natural frequencies, maximum stress at the root, stage weight and tip 
displacements. Geometric and material perturbations have been conducted for 
an SSME blade using a special purpose code based on finite element analysis. 

The geometric perturbations which simulate the natural perturbations under 
operating and/or fabrication conditions, are generated by randomly perturbing 
the x, y, and z coordinates of all the nodes of the finite element mesh. The 
material perturbations were generated by randomly perturbing the material pro- 
perties of all finite elements. These perturbations represent In part material 
variation resulting from the fabrication process, and/or any other local 
Irregularities. 
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Probabilistic models have been developed to predict the structural response 
by perturbatlng means and variances. Since means are not significant, prob- 
abilistic models based only on the significant variances have also been devel- 
oped. In addition, probability distributions for the structural response have 
been developed. These distributions provide an assessment of variation In the 
structural response for selected geometric and material properties perturba- 
tions. These provide an estimate of the probability of getting a certain res- 
ponse for a given Input. Some representative results are shown In figure 14. 
Statistical tests and methods were used to check the developed models. These 
tests Indicated that the developed models are good fits. T-tests were used to 
Identify the significant variances of perturbations. F-tests, plots, and 
autocorrelations of residuals were used to check the goodness-of-f It tests. 


SUMMARY 

NASA Lewis Research Center Is currently developing probabilistic structural 
analysis methods for select SSME structural components. Briefly, the develop- 
ment consists of the following program elements: (1) composite load spectra, 

(2) probabilistic structural analysis methods, (3) probabilistic finite element 
theory - new variational principles, and (4) probabilistic structural analysis 
application. The development of the probabilistic structural analysis method- 
ology Is a joint effort of NASA Lewis In-house research, contract, and grants. 
The research activities have led to significant technical progress In several 
Important aspects of probabilistic structural analysis. The significant tech- 
nical accomplishments to date demonstrate that structural analyses can be for- 
mulated using probabilistic methods where all the participating structural- 
component-descriptor parameters and variables, and loading conditions are 
defined probabilistically. An early version of a structural analysis computer 
code (NESSUS) based on probabilistic finite elements has been completed and Is 
currently used to analyze high pressure turbopump blades. Also the Individual 
loads and composite load spectra can be probabilistically simulated using 
probabilistic methods with progressive levels of sophistication, limited 
available data, expert opinion, and an expert system driven computer code. 
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TABLE 1. - INDIVIDUAL LOAD SUMMARY 


Individual loads 

How determined 

Form available 

Degree of 
certainty 

How used 

Centrifugal 

Measured 

Duty cycle 
and limits 

High 

1 

Input to structural analysis 
as values, limit cases or 
duty cycle 

Static pressure 
and aP 

Predicted from 
turbine analysis 

Pressure profiles 
at streamlines on 
air foil 9 

Total magnitude 
high (torque) dis- 
tribution moderate 

Input to structural analysis 
as limit cases or duty cycle 

Dynamic pressure 

Predicted from 
turbine analysis 
and cascade 
estimates or 
tests 

Pressure profile 
scaled to forcing 
function shape 9 

Moderate to low 

Input to forced vibration 
analysis, Campbell diagram 
1 imits 

Temperature 

Predicted from 
turbine plus 
thermal analysis 

Temperatures at 
location throughout 
blade (computerized 
dataset) 9 

Steady state-high 
transient-low 

Input to structural analysis 

Debris 

Predicted, based 
on history of 
incidents 

Particle size and 
velocity 

Low 

Sensitivity analysis or 
incident investigation 


a Low frequency and transient 


TABLE 2. - SUMMARY MATRIX OF INDIVIDUAL LOAD VERSUS COMPONENT 


Individual load 

Turbine 
bl ade 

Transfer 

duct 

Lox 

post 

HPOTPDD 

Load form 

Static pressure 

X 

X 

X 

X 

Duty cycle 9 

Dynamic pressure 






Chugging (transient) 

- 

X 

- 

- 

AMS. statos 

Tubulence 

Sinusidal 






(repeated pulse) 

X 

X 

- 

- 

AMS. PSD. statos 

Random 

- 

X 

X 

X 

AMS. PSD. 

Centrifugal 

X 

- 

- 

- 

Duty cycle 9 

Temperature 

X 

X 

X 

X 

Duty cycle 9 

Structural vibration 
Transient 






Sideload 

- 

X 

X 

_ 

AMS. statos 

Pops 

- 

X 

X 

- 

AMS. statos 

Steady state 






Sine 

- 

X 

X 

X 

AMS. PSD. statos 

Random 

- 

X 

X 

X 

AMS. statos 

Debris 

X 

X 

X 

- 

History 

Rubbing 

X 

- 

- 

- 

Expert opinion 

Instal lation 

- 

- 

X 

X 

Expert opinion 

Fab 

X 

X 

X 

X 


Friction 

X 

X 

X 

_ 

PSEUDO load 

Tolerances 

X 

X 

X 

X 



a Low frequency and transient 
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Figure 1 . - PSAM will be initially developed for select SSME structural components. 
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Figure 3. - Composite load spectra flow chart/task features summary. 
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Figure 4. - High-pressure fuel turbine transient characteristics. 
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Figure 7. - Key technical features of composite load code. 
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Figure 8. - Composite load spectra simulation using expert systems. 



Figure 9. - PSAM integrates existing technologies. 
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PFEM EXAMPLE PROBLEM - CURVED SHELL 
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(NODE 1) 
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COEFFICIENT OF VON MISES STRESS AT NODE 1 


RANDOM VARIABLE MEAN VARIATION (COV) DISTRIBUTION INDEPENDENT 


SHELL THICKNESS 

0.15 IN. 

0.05 
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YES 

FACTOR 





TEMPERATURE 
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NO 

PRESSURE 
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MATERIAL YIELD 
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STRENGTH 







Figure 12. - Probabilistic finite element stress analysis. 
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Figure 13. - Variational approach to probabilistic finite elements. 
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Figure 19. - Probabilistic structural analysis can be used to explain SSME turbopump blade freouency distributions. 
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